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Abstract 

A  let  of  general  sequential  filter  equations  is 
derived  for  nonlinear  system  dynamics  and  a  non¬ 
linear  observation  model,  but  is  obtained  with  the 
assumption  of  a  linear  estimator.  These  equa¬ 
tions,  which  are  based  upon  previously  developed 
formulas,  include  the  colored  measurement  noise 
statistics  and  the  statistics  of  nonestimated  model 
parameter  errors.  Several  simulations  made  with 
this  filter  are  compared  with  the  simulation 
results  of  an  estimation  procedure  constructed 
with  the  utilization  of  the  standard  white  noise 
assumptions.  The  difference  between  the  white 
noise  filter  results  and  the  colored  noise  filter 
results  is  found  to  be  minimal.  Instability  occur¬ 
red  when  the  statistics  of  the  effective  exhaust- 
velocity  in  the  acceleration  model  were  not  proper¬ 
ly  accounted  for. 
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L  Introduction 


During  radio-guided  launch  miaeiona,  the 
angle  meaaurement  noiae  of  the  guidance  radar 
becomea  increaalngly  botheraome  aa  the  radar 
line-of-aight  tenda  toward  low  elevation  anglea. 
Significantly,  a  greater  portion  of  the  launch 
guidance  ha  a  been  occurring  at  low  elevation  anglea 
aa  a  conaequence  of  the  more  intenaive  uae  of 
larger  booat  vehiclea,  which  have  become  avai¬ 
lable  to  place  heavier  payloada  in  orbit.  The 
botheraome  component  of  the  angle  meaaurement 
noiae  ia  due  to  random  fluctuationa  in  the  tropo- 
apheric  index  of  refraction.  Theae  fluctuationa 
are  moat  noticeable  in  the  denaer  portion  of  the 
atmosphere  and  are  directly  attributable  to  the 
random  nature  of  temperature,  preaaure,  and 
humidity  variations.  These  variations  are  slower 
and  more  severe  in  the  low-altitude  regions  and 
thus  affect  not  only  the  magnitude  but  also  the 
autocorrelation  of  the  angle  errors  when  the  radar 
energy  traverses  these  regions.  In  fact,  angle 
noiae  correlation  times  of  several  seconds  are 
common  at  radar  elevation  anglea  of  a  few  degrees. 
In  addition,  the  present  trend  in  guidance  filter¬ 
ing  has  been  toward  more  intensive  use  of  Kalman 
filtering  techniques.  Since  filters  of  this  type 
require  that  the  measurement  noiae  characteristics 
be  modeled,  it  is  appropriate  to  investigate  the 
behavior  of  such  filters  when  various  assumptions 
are  made  concerning  the  noise  correlations. 

Error  sources  other  than  measurement  noise 
affect  the  guidance  filter  behavior  and  are  of 
prime  concern  in  a  radio-guided  launch  mission. 

A  filter  constructed  to  minimise  the  effects  of 
meaaurement  noise  will  itself  introduce  errors 
into  the  estimate  that  arise  from  inaccuracies  in 
specifying  the  mathematical  model  of  the  system. 
These  modeling  errors  reflect  the  limits  inhuman 
ability  to  describe  the  real  world  in  exact  mathe¬ 
matical  terms.  Errors  of  this  sort  can  have  an 
accumulative  effect  over  a  period  of  time  and  can 
cause  the  estimate  to  diverge.  It  is  assumed  for 
the  purposes  of  this  paper  that  the  mathematical 
form  of  the  system  model  is  known  or  can  be 
adequately  represented  by  some  series  expansion; 
however,  the  values  of  certain  parameters  of  the 
model  are  in  error.  This  assumption  differs 
significantly  from  the  situation  where  the  mathe¬ 
matical  form  itself  is  uncertain.  Examples  of 
model  parameter  uncertainties  are  survey  uncer¬ 
tainties,  measurement  bias  errors,  and  para¬ 
meter  errors  such  as  the  uncertainty  in  the 
numerical  value  of  the  gravitational  constant. 


In  this  paper  a  derivation  is  given  of  a  general 
aequential  filter  which  utilizes  the  differencing 
scheme  of  Refs.  1  and  2  to  account  for  sequentially 
correlated  noise.  Nonestimated  model  parameter 
uncertainties  in  both  the  state  dynamic  model  and 
the  measurement  model  are  accounted  for  ia  the 
derivation. 

The  various  features  of  the  general  filter 
were  programmed  in  a  simulation  of  a  launch 
guidance  mission  and  the  results  are  given  here. 
The  primary  purpose  of  the  simulation  was  to 
examine  the  behavior  of  the  general  filter  and  to 
determine  whether  its  performance  was  such  as  to 
warrant  an  increase  in  the  complexity  of  existing 
guidance  programs. 

II.  The  Optimum  Filter 

Before  proceeding  into  the  mathematical 
details  of  the  derivation,  a  qualitative  description 
of  the  filter  will  be  given  to  help  in  visualising  the 
filtering  process  in  the  presence  of  correlated 
measurement  noise  and  model  parameter  errors. 
An  optimum  linear  sequential  filter  operates  on 
data  samples  that  are  available  at  discrete  values 
of  time.  Its  output  at  a  particular  time  is  an  esti¬ 
mate  of  the  system  state  based  on  all  information 
available  up  to  that  time.  A  block  diagram  of  a 
sequential  filter  is  shown  in  Fig.  1.  It  is  optimum 
in  a  least  squares  or  minimum  variance  sense, 
providing  that  the  Gaussian  and  linear  propagation 
properties  of  the  errors  hold;  and  it  is  linear 
because  the  estimates  are  formed  from  a  linear 
combination  of  the  observation  data. 


At  the  beginning  of  each  computation  cycle 
there  are  available  an  estimate  of  the  state  vector 
and  a  covariance  matrix  of  errors  in  the  estimate. 
Also  input  at  this  point  are  a  set  of  measurements 
and  a  covariance  matrix  of  measurement  errors. 
These  data  -•  toe  state  estimate  and  the  measure¬ 
ment  set  --  are  weighted  according  to  the  relative 
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magnitude  of  their  respective  covariances,  data 
with  larger  covariances  being  weighted  less  than 
those  with  smaller  covariances,  and  are  appro¬ 
priately  transformed  to  form  a  new  state  vector 
estimate.  Next,  a  new  state  vector  estimate  error 
covariance  matrix  is  computed  by  propagating  the 
previous  state  vector  estimate  error  covariance 
matrix  and  the  measurement  error  covariance 
matrix  through  the  filtering  equations.  Finally, 
the  state  vector  estimate  and  the  state  vector 
estimate  error  covariance  matrix  are  propagated 
to  the  beginning  of  the  next  computation  cycle  by 
the  filter  dynamic  and  error  propagation  models. 

Consider  now  the  situation  where  the  mea¬ 
surement  noise  is  highly  correlated,  as  would  be 
the  case  in  low-elevation  radar  tracking.  The 
estimate  errors  will  be  correlated  with  the 
measurement  noise  through  previous  measure¬ 
ments.  A  filter  based  on  white  noise  assumptions 
will  ignore  this  fact  and  improperly  propagate  the 
covariance  matrices,  consequently  causing  im¬ 
proper  weighting  of  the  data.  For  example,  sup¬ 
pose  that  the  range  from  a  radar  to  a  fixed  point 
were  to  be  determined  from  range  measurement 
samples  and  that  these  measurements  possessed 
noise  with  long-term  drift  characteristics.  A 
white  noise  filter  would  treat  the  measurement 
noise  as  being  statistically  independent  from 
sample  to  sample  with  a  statistical  mean  that 
tended  toward  aero  as  the  number  of  samples 
increased,  when  in  actuality  the  measurement 
error  would  behave  more  as  a  bias  if  the  total 
time  span  under  consideration  is  short. 

One  means  of  accounting  for  the  presence  of 
colored  noise  is  to  augment  the  state  vector  to 
include  the  states  of  a  hypothetical  shaping  filter 
whose  input  is  white  noise  and  whose  output  is 
noise  of  the  required  correlations.  This  approach 
is  often  prohibitive  primarily  because  of  computer 
storage  limitations,  since  for  each  increase  in  the 
state  vector  dimensions  there  is  a  corresponding 
increase  in  the  dimensions  of  the  state  covariance 
matrices  and  in  the  complexity  of  the  ensuing 
matrix  arithmetic.  Also,  the  state  augmentation 
scheme  can  result  in  ill-conditioned  matrices  (1). 

A  second  method  which  utilises  a  measurement 
differencing  technique^)*  ^  reduces  the  imple¬ 
mentations!  difficulties  associated  with  state  aug¬ 
mentation  and  at  the  same  time  alleviate  the 
threat  of  matrix  inversion  problems.  Because  of 
these  advantages,  the  second  method  was  used  in 
the  launch  guidance  simulation.  Recently,  a  third 
method  was  developed  which  does  not  depend  upon 
state  augmentation  or  measurement  differencing. 

It  was,  however,  developed  too  late  to  obtain 
results  for  the  paper,  and  hence  will  be  discussed 
in  a  subsequent  report^ ). 

Basically,  there  are  three  ways  in  which  the 
model  parameter  errors  can  be  accounted  for  ui 
the  estimation  process.  First,  they  can  be  ig¬ 
nored  completely.  In  this  case  the  state  vector 
estimate  error  covariance  matrix,  since  it  will 
not  show  the  effect  of  propagating  the  parameter 
errors,  will  be  smaller  than  it  should  be.  Conse¬ 
quently,  the  state  vector  estimate  that  is  present 
at  the  beginning  of  each  computation  cycle  will  be 
weighted  more  and  the  measurements  weighted 
less  than  if  the  effects  of  the  parameter  errors 
were  accounted  for. 

Secondly,  the  parameters  can  be  estimated 
along  with  the  rest  of  the  state  vector  elements  in 
hopes  of  driving  the  bothersome  errors  to  smaller 


values.  In  any  event,  the  parameter  error  co- 
variances  will  be  included  in  the  state  vector 
estimate  error  covariance  matrix  and  will  hence 
Influence  the  computation  of  the  filter  weights. 

That  is,  even  if  the  parameter  uncertainties  are 
not  reduced,  they  will  be  accounted  for  in  the 
filter  error  analyais. 

A  third  way  to  account  for  the  model  para¬ 
meter  uncertainties  is  similar  in  some  respects 
to  the  second  method.  In  this  method  the  error 
covariances  of  the  parameter  errors  are  included 
in  the  filter  covariance  matrix  propagations!4), 
but  unlike  the  second  method  the  parameters  are 
not  solved  for,  and  thus  the  parameter  errors 
never  decrease.  The  procedure  is  equivalent  to 
the  situation  where  the  parameter  is  solved  for 
as  in  the  second  method  but  where  poor  parameter 
estimates  are  obtained.  Attempts  have  been  made 
to  approximate  this  third  approach  by  the  addition 
of  an  artificial  "state  noise"  covariance  matrix  to 
the  state  error  covariance  matrix  once  during 
each  computation  cycle.  Although  such  an  ap¬ 
proach  is  not  a  systematic  one,  it  is  obviously 
better  than  ignoring  the  parameter  errors  as  the 
first  method  does  since  it  prevents  the  state  error 
covariance  matrix  from  becoming  too  small. 
Nevertheless,  the  correlation  between  the  state 
estimate  error  and  the  parameter  error  is  ignored 
inasmuch  as  the  state  noise  approach  is  based  on 
the  assumption  of  white  noise  parameter  errors. 

IIL  Derivation  of  the  General  Filter 

Let  the  state  dynamics  be  described  by 

xn  =  <P(xn.1(  P.Un)  (1) 

and  the  relation  of  the  measurement  to  the  states  by 
zn  =  h(Xn,T)  +  Nn  (2) 

Let  the  measurement  noise  be  generated  by  the 
Markov  process 

Nn  =  ’Vl  +*Wn  <3> 

The  filter  dynamic  model  is  assumed  to  be  of 
the  form 

*n/n-l  =  ^n-l/n-1’  <4> 

and  the  filter  measurement  model  of  the  form 

*'n  *  h^n/n-ri)  (5) 

Let  the  measurements  be  differenced  to  give 
an  effective  measurement 


C  =  s  -  p  z  , 
n  n  n-1 

The  measurement  estimate  is 

{  *  ^n/n-l'^-^n-l/n-l’T) 

■  h.(in-l/n-l’*> 


(6) 


(7) 
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Assume  that  the  state  vector  estimates  are 
obtained  from  a  linear  filter  of  the  form 

*n/n  =  *n/n-1+lSA-<n1  <8> 

The  error  in  the  estimate  is 

8*n/„  "  Xn  '  *n/n  =  6^n/n-l  ‘  K„<  <9> 


Eft=h8^  -  HaP„-lHaT  +  ^  +  HnrarT«n 


+  H  DVT  ,HT  +  H  V  DTHT 
n  n-1  a  a  n-1  n 


+  H  DSDTHT  +  alJ  .  HT  +H  L  .  AT 
n  n  n-1  a  a  n-1 


8*n/n-l  "  Xn  *  *»/n-l  <10> 

8^a£n-^n  <U> 

The  filter  described  by  the  above  equations  is 
shown  in  the  block  diagram  of  Fig.  1. 

By  the  Wiener-Hopf  Equation  (see  Appendix) 
the  filter  gain  is  obtained  as 

Kn'H/n-l60|«8^lTf1  <12> 

To  obtain  an  explicit  expression  for  K  it  is  first 
necessary  to  expand  Eq.  (10).  Thus  n 

8*n/n-l  =’»Xn-l'^n|-<1/„-l^1  (13) 


P  =  (*-K  H  )P  ,(*-K  H  )T  +  K  qRqTKT 
n  n  a  n-l  n  a  n  ’  n 


+  (I-K  H  )rOTT(l-K  H  )T 

n  n  n  n 


+  (K-K  H  )DSDT(I-K  H  )T 
n  n  n  n 

+  (♦  -  K  H  ) V  .DT(I-KH  )T 
n  a  n-l  n  n 

+  KnHn  ®  Vn  - 1  ”KnHa  ^ T 

+  K  ABATKT  -  (*-K  H  )L  ATkJ 
n  n  n  a  n-i  n 

-  K  AlJ  (t-K  H  )T 
n  n- i  n  a 


which,  if  one  assumes  that  the  estimate  is  within 
a  linear  expansion  of  the  state,  becomes 


L  =  (t-K  H  )L  -  K  AB 
n  n  a  n-1  n 


Wl  =48Xn-l/„-l+D8P+I^n 


V  =  (t-K  H  W  +  (I-K  H JDS  (20) 

n  n  a  n-l  n  n 


Similarly,  Eq.  (11)  can  be  written  as 

it  =  (H  t-pH  ,)SX  .  +  H  Dtp 

n  '  n  n-1'  n-1 /n-1  n  r 


Equations  (16)  through  (Z0)  along  with  Eq.  (1Z) 
constitute  the  general  sequential  filter. 

IV.  Application  of  the  Filter 


+  AZT  +  H  r u  +  qW 
n  n  n  n 


=  H  6X  . +H  D8p  +  AZT  +  H  I\i  +  qW 

a  n-l/n-1  n  r  n  n  ^  n 


After  one  uses  Eqs.  (14)  and  (IS)  in  Eq.  (12),  and 
performs  the  indicated  expected  values,  K  can 
be  determined  from  n 

E[«X  .  1«£T1  =  tP  HT  +  t  V  ,DTHT 

n/n-1  n  n-1  a  n-1  n 

+  4L  ,AT  +  DVT  ,HT 
n-1  n-1  a 


+  dsdtht  +  rorT  ht 

n  n 


In  radio-guided  ascent  missions  a  tracking 
radar  measures  the  position  and/or  velocity  of  the 
launch  vehicle  (see  Fig.  Z).  These  measurements 
are  processed  by  a  guidance  computer  to  obtain 
an  estimate  of  the  vehicle  state.  The  computer 
determines,  according  to  a  programmed  guidance 
law,  what  vehicle  steering  is  necessary  to  satis¬ 
fy  the  mission  requirements;  and  steering  com¬ 
mands  are  transmitted  to  the  vehicle  via  the 
guidance  radar  link.  In  addition,  the  radar  trans¬ 
mits  discrete  commands  such  as  engine  cutoff  as 
soon  as  the  computer  determines  the  necessity  of 
such  events. 

A  simulation  of  a  typical  launch  guided  mis¬ 
sion  was  constructed,  as  shown  schematically  in 
Fig.  3,  in  which  a  vehicle  was  guided  during  its 
last  two  stages.  The  purpose  of  the  simulation 
was  to  assess  the  performance  of  the  general  se¬ 
quential  filter  in  a  simulated  radar  noise  environ¬ 
ment,  and  in  the  presence  of  uncertainties  in  both 


Fig.  2.  Launch  Guidance  Radar  System 

the  filter  state  model  parameters  and  filter  mea¬ 
surement  parameters.  The  vehicle  was  assumed 
to  be  guided  perfectly,  hence  the  state  vector  in 
the  simulation  was  independent  of  the  filter  esti¬ 
mates.  The  performance  factors  were  the  com¬ 
puted  residuals  of  the  state  estimation  error  and 
the  computed  state  estimation  error  covariance 
matrix.  A  position  radar  measuring  range,  azi¬ 
muth,  and  elevation  was  simulated.  The  simula¬ 
tion  error  model  added  random  noise  to  the  simu¬ 
lated  measurements  which  possessed  the  following 
one-half  second  correlations: 

range  noise  correlation  =  0.  1 
azimuth  noise  correlation  =  0.  95 
elevation  noise  correlation  =  0. 97 

The  estimated  state  was  a  12-element  vector  con¬ 
sisting  of  the  following  quantities: 

x,  y,  s  -  a  3-element  ECI  Cartesian 
position  vector 

x.  y,  i  -  a  3-element  ECI  Cartesian 
velocity  vector 

c*  -  effective  engine  exhaust  velocity 
tm  -  engine  mass  depletion  time 

i  -  x-component  of  a  roll  axis  oriented 
unit  vector 

i  -  z- component  of  a  roll  axis  oriented 
unit  vector 

flx  -  pitch  axis  drift  rate 
'  y*w  axis  drift  rate 

The  following  7  simulations  were  performed: 

s.  The  filter  error  model  was  based  on  an 
assumption  of  white  measurement  noise 
(white  noise  filter).  The  simulation 
error  model  generated  measurement 
noise  with  the  correlations  given  above. 


Fig.  3.  Simulation  Block  Diagram 

b.  The  filter  error  model  was  based  on  an 
assumption  of  correlated  measurement 
noise  with  the  correlations  given  above 
(colored  noise  filter).  The  simulation 
error  model  generated  measurement 
noise  with  these  same  correlations. 

c.  The  filter  was  a  white  noise  filter,  and  the 
simulation  error  model  generated  bias 
errors  onlv. 

d.  The  filter  error  model  was  based  on  an 
assumption  of  measurement  bias  errors, 
and  the  simulation  error  model  generated 
bias  errors  only. 

e.  The  filter  was  a  white  noise  filter,  the 
simulated  radar  measurements  were 
error-free,  and  c*  was  omitted  from  the 
estimated  state  vector.  No  compensa¬ 
tion  for  c*  uncertainties  was  attempted. 

f.  The  filter  was  a  white  noise  filter  and  the 
simulated  radar  measurements  were 
error-free,  c*  was  included  in  the  esti¬ 
mation  vector. 

g.  This  simulation  was  identical  with  simu¬ 
lation  e  except  that  the  c*  error  statis¬ 
tics  were  included  in  the  filter  equations. 

V.  Results 

Typical  results  of  the  first  simulation  (the 
white  noise  filter)  are  given  in  Fig.  4a  and  can  be 
compared  with  the  results  of  the  second  simula¬ 
tion  (the  colored  noise  filter)  given  in  Fig.  4b. 

Goth  figures  show  the  x-components  of  the  velocity 
estimate  residuals.  In  the  same  plots  the  x- 
velocity  standard  deviation  envelopes  obtained 
from  the  state  estimate  error  covariance  matrix 
of  the  filter  are  also  shown.  Comparison  of  the 
white  noise  filter  velocity  residuals  of  Fi?.  4a 
with  the  velocity  residuals  of  the  colored  noise 
filter  in  Fig.  4b  show  that,  although  there  is  an 
improvement  in  the  estimate,  the  improvement  is 
small,  because  the  average  deviation  of  the  resi¬ 
dual  about  the  zero  mean  value  is  about  the  same 
in  both  simulations.  There  is,  for  example,  only 
about  a  2  ft/sec  difference  in  the  velocity  residual 
at  burnout.  The  filter  value  of  the  estimate  error 
standard  deviations  as  shown  by  the  standard 
deviation  envelopes,  however,  is  larger  in  the 
colored  noise  case  and  is  more  representative  of 
the  estimate  errors. 
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Fig.  4a.  Velocity  Residual,  White 
Noise  Filter 


Fig.  5a.  Velocity  Residual,  White 
Noise  Filter 
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Fig.  4b.  Velocity  Residual,  Colored 
Noise  Filter 


Fig.  Sb.  Velocity  Residual, 
Bias  Filter 


The  results  of  simulations  c  and  d  showing 
the  effects  of  measurement  bias  errors  on  the 
output  of  the  white  noise  filter  and  the  bias  filter 
are  given  in  Figs.  5a  and  5b.  A  small  difference 
can  be  seen  between  the  white  noise  filter  results 
and  the  bias  filter  results.  The  computed  stand¬ 
ard  deviations  were  larger  when  the  bias  filter 
was  used  as  illustrated  by  the  larger  values  of  the 
envelope. 


Figures  6  and  7  show  the  results  of  the  simu¬ 
lation  e  in  which  c*  was  omitted  from  the  state 
estimation  vector  and  the  c*  uncertainties  were 
ignored.  Figure  6  shows  that  the  filter  velocity 
estimates  become  unstable  after  only  about  60  sec. 
The  filter  gains  plotted  in  Fig.  7  are  the  gains 
relating  the  position  estimates  to  the  transformed 
radar  measurements;  the  radar  measurements 
are  transformed  into  the  same  coordinates  as  the 
state  vector  position  components.  The  gains  in 
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Fig.  7  become  ieee  then  0.4  after  the  first  10  sec 
of  filtering  which  implies  that  the  a  priori  esti¬ 
mate  —  the  state  estimate  present  at  the  beginning 
of  each  computation  cycle  --  receives  a  greater 
share  of  the  weighting  than  do  the  measurements. 
The  somewhat  wild  behavior  of  the  gains  during 
the  second  stage  is  attributed  to  the  unstable 
behavior  of  the  estimates  that  are  used  to  evaluate 
the  partial  derivatives  which,  in  turn,  affect  the 
gain  computations. 


Fig.  o.  Velocity  Residual;  White  Noise  Filter, 
c*  not  Estimated 


Fig.  7.  Filter  Gain,  White  Noise  Filter, 
c*  not  Solved  for,  c*  Errors 
Ignored 


results  of  simulation  g  in  which  c*  was  not  esti¬ 
mated  but  the  c*  errors  accounted  for.  The  posi¬ 
tion,  velocity,  and  attitude  estimate  residuals  are 
nearly  the  same  in  both  simulations.  Note  that  the 
filter  position  gains  of  Fig.  9b  are  always  greater 
than  those  of  Fig.  7  during  the  first  stage  demon¬ 
strating  that  the  inclusion  of  the  c *  statistics  in 
the  filter  results  in  greater  weighting  of  the  mea¬ 
surements  than  when  they  are  ignored.  The  gains 
of  Fig.  9b  also  lie  above  those  in  Fig.  9a.  Since 
the  c*  estimate  error  decreased  with  time  in  the 
simulation  whose  results  are  shown  in  Fig.  9b, 
they  therefore  had  a  smaller  effect  when  the  state 
estimate  was  propagated.  This  error  decrease  is 
implicitly  accounted  for  in  the  filter  equations  in 
the  covariance  matrix  computations. 


Fig.  8a.  Velocity  Residual,  White 
Noise  Filter 


Figure  8a  shows  the  results  of  simulation  f  in 
which  c*  was  estimated,  and  Fig.  8b  shows  the 


Fig.  8b.  Velocity  Residual,  c*  Errors 
Accounted  for 
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Fig.  9a.  Filter  Gain,  White  Noise 
Filter,  c*  Solved  for 


Fig.  9b.  Filter  Gain,  White  Noise  Filter, 
c*  not  Solved  for,  c*  Errors 
Accounted  for 


VL  Conclusions 

The  simulation  results  indicate,  as  would  be 
expected,  that  the  colored  noise  filter  perfor¬ 
mance  in  a  colored  noise  environment  is  superior 
to  the  performance  of  a  white  noise  filter  in  the 
same  noise  environment  since  the  colored  noise 
filter  estimation  errors  are  smaller.  However, 
the  degree  of  improvement  is  minimal,  and  in  a 
practical  sense,  the  improvement  would  not  war¬ 
rant  the  increased  complexity  of  the  guidance 
equations.  Similarly,  the  bias  filter  does  not, 
for  this  application,  exhibit  enough  performance 
improvement  to  justify  its  implementation. 


State  vector  estimates  which  do  not  contain  c* 
as  one  of  its  elements  can  be  satisfactorily  ob¬ 
tained  if  the  c*  error  statistics  are  properly 
accounted  for  in  the  filter  equations.  Some  form 
of  compensation  for  the  omission  of  this  parameter 
from  the  estimate  is  mandatory  if  the  filter  is  to 
remain  stable. 

In  general,  the  results  demonstrate  how  cau¬ 
tion  must  accompany  the  desire  for  optimality  in 
a  practical  application  of  theoretical  results.  It 
also  demonstrates  how  the  optimal  approach,  whan 
used  as  an  analysis  tool,  provides  useful  goals  or 
limitations  upon  which  a  practical  design  can  be 
based. 

Appendix 

The  Wiener-Hopf  Equation 

Define  the  inner  product  of  two  real  valued 
vectors  a  and  b  to  be 

<a,  b>  ■  XX  E  (abT)  X2  (A-l) 

where  a  is  an  n-vector  and  b  is  an  m-vector.  Let 
X.  be  any  real  n-vector  and  X,  any  real  m-vector 
stlch  that  2 

X  j  C  X^  for  n  <  m 

Xj  0  X2  for  m  <  n 

X|  =  Xj  for  m  =  n 

An  inner  product  of  two  real  vectors  must 
satisfy 

a.  <orX  +  By,  a  >  =  or<  X,  s  >  +  B<  y,  s> 
where  a  and  B  are  real  scalars. 

b.  <X,  X>2  0;  <  X,  X>  =  0  ^X  =  0 

Satisfaction  of  condition  a  by  Eq.  (A-l)  is  easily 
demonstrated.  Condition  b  merely  states,  for  the 
inner  product  defined,  the  positive -definite  pro¬ 
perty  of  a  covariance  matrix. 

A 

The  vector  5X  ,  of  Eq.  (9 )  in  the  text  is  the 
estimation  error  aTra  represents  the  additional 
information  required  to  specify  the  state  exactly. 
The  vector  it  represents  the  information,  in 
addition  to  that  obtained  from  past  data,  provided 
by  the  most  recent  observation.  If  iXn/n  is  an 
optimal  estimate  in  a  least  squares  sense, 
is  orthogonal  to  it  ;  i.  e., 

<4*n/n-4*h>=  0  (A-2) 

But  if 

xj  E(abX)  X2  =  0 

for  all  X.  and  X, ,  then 
•  1  Z 

E  (abT)  =  0 
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•o  that 


»♦  *  T  -  t 


*  0 

*  *»*«/«-!  •  -  0 

If  on*  »olv«»  (or  K,  th«  Wi«n«r-Hop(  Equation  U 

obtained  aa 


B  •  Etsia-P] 

.  x^£  ,  , 

D  =  |S|  *n/n*1 

dp'p  =  p 


E  =  the  expectation  operator 


Vn  "  *»*./« 


W  a  a  white  noiae  vector 

n 

X_  ■  the  atate  vector 

n 

X  /i  a  the  eatlmate  of  the  atate  vector  at  tR 
'  given  all  meaaurementa  up  to  and 
including  that  at  tj 

*  a  the  meaaurement  vector 


=  effective  meaaurement  after  differencing 

=  estimated  effective  meaaurement  computed 
£rom  n/n-1 

U  =  white  proceaa  noise  vector 
n 

X  =  4 

r  ■  |S) 

ap  P  =  P 


"a  =  «h*-'*«n.l 

A 

H  -^-lX=an/n*1 

"n  *  5T3C  1  T  =  T 


Kr  a  the  filter  gain  matrix 

Ln*  *4*n/n4*T1 

n  =  subscript  associating  variable  with  the 
time  point 

n/n-1  =  subscript  indicating  that  the  variable  ia 
evaluated  at  the  n*“  time  point  and  is 
based  upon  measurements  up  to  and 
including  the  (n-1)** 


N  a  vector  of  measurement  noise 

n 

p  a  dynamic  model  parameter  vector 

p  a  an  estimate  of  p 

*  * 

«p  =  p  -  p 

Q  ■  E&i  ujl 

n  d 

R  -  E[W  wj] 

n  n 

S  ■  Eftp  4p^l 

T  a  measurement  model  parameter  vector 


p  =  correlation  coefficient  matrix  representing  the 


correlation  between  N  and  N  . 

n  n-i 


1  P  -  P 
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